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Abstract. Stark deceleration has been utilized for slowing and trapping several species of neutral, ground- 
state polar molecules generated in a supersonic beam expansion. Due to the finite physical dimension 
of the electrode array and practical limitations of the applicable electric fields, only molecules within a 
specific range of velocities and positions can be efficiently slowed and trapped. These constraints result 
in a restricted phase space acceptance of the decelerator in directions both transverse and parallel to the 
molecular beam axis; hence, careful modeling is required for understanding and achieving efficient Stark 
decelerator operation. We present work on slowing of the hydroxyl radical (OH) elucidating the physics 
controlling the evolution of the molecular phase space packets both with experimental results and model 
calculations. From these results we deduce experimental conditions necessary for efficient operation of a 
Stark decelerator. 

PACS. 32.60. +i Zeeman and Stark effects - 39.10+j Atomic and molecular beam sources and techniques 



1 Introduction 

The study of ultra-cold polar molecules is a rapidly emerg- 
ing field as evidenced by this special journal edition, and 
the large number of new experiments being undertaken 
PQEBSEISlini • The permanent electric dipole moment, /i, 
possessed by polar molecules allows for a new type of inter- 
action in the ultra-cold environment. This electric dipole- 
dipole interaction (and control over it) should give rise to 
unique physics and chemistry including novel cold collision 
dynamics quantum information processing jHj, and a 
second order phase transition analogous to spontaneous 
magnetization of a ferromagnet. When the translational 
energy of colliding polar molecules becomes comparable 
to their dipole-dipole interaction energy, the molecules 
can dramatically influence each other's trajectory. Lack 
of spherical symmetry in the interaction causes collid- 
ing molecules to be attracted or repelled depending on 
their relative orientation. Thus, an external electric field, 
which orients the molecules, may have a profound effect 
on the molecular interactions leading to large changes in 
observed scattering rates [Zj . Another interesting possibil- 
ity is the observation of criticality in a two or three (3D) 
dimensional lattice of dipolar molecules. For a typical po- 
lar molecule (~ 2 Debye) in a 3D optical lattice with 1 
m spacing, the critical temperature for the second order 
phase transition (spontaneous polarization) is estimated 
by a 3D Ising model as T c - 3/i 2 / (27re k B R 3 ) = 200 nK, 
where R is the dipole-dipole separation, eo is the permit- 
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tivity of free space, and fc# is Boltzmann's constant jHj. By 
using extremely polar molecules such as those predicted 
to occur from photo-association of free radicals with Ru- 
bidium ( \i ~ 10 Debye) [TOl, or by decreasing the dipole- 
dipole separation through simply trapping polar molecules 
at a sufficient density, the critical temperature may be 
raised by at least one order of magnitude. It is interesting 
to note that at temperatures significantly higher than this 
regime an analog to adiabatic demagnetization refrigera- 
tion [ll] may become possible, leading to a new method 
for further cooling polar molecules, so-called "paraelectric 
cooling" H2|. 



The technique of Stark deceleration is especially well- 
suited for production of cold polar molecules since it uti- 
lizes the Stark shift associated with polar molecules to 
remove the mean translational energy from a supersonic 
beam of molecules. Supersonic molecular beams, utilized 
extensively in physical chemistry, are capable of producing 
intense pulses of molecules with mean velocities of a few 
hundred meters per second and a small velocity spread 
about this mean. Therefore, in the moving frame of the 
molecular pulse the associated velocity spread corresponds 
to a low temperature distribution IK). Stark decelera- 
tion conservatively removes this mean velocity, leaving be- 
hind a cold molecule distribution in the laboratory frame. 
To date this technique has been employed to produce slow, 
cold beams of CO 0, OH jSJ, ND3 resulting in a 
trapped sample of 10 4 molecules at 25 mK in the case 
of ND3. Slowing of YbF in an alternate-gradient deceler- 
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ator has also been demonstrated in a proof of principle 
experiment jHj. 

Due to practical limitations of both the physical di- 
mensions of a Stark decelerator and the applicable electric 
fields, only molecules within a specific range of velocities 
and positions can be efficiently slowed and trapped. These 
constraints result in a restricted phase space acceptance of 
the decelerator in directions both transverse and parallel 
to the molecular beam axis. Thus, for efficient operation 
care must be taken to match the phase space distribution 
of the supersonic beam source (emittance) to the phase 
space acceptance of the Stark decelerator. This article de- 
tails work from both experiment and model, describing the 
process of phase space matching between the supersonic 
beam source and the molecular decelerator. The focus of 
the current article is to present a clear understanding of 
the decelerator efficiency and describe explicitly the exper- 
imental requirements of maximizing it. For experimenter 
new to this field, this article will serve as a useful and 
practical guide for design of future Stark deceleration ex- 
periments. 

Our experiment centers on the deceleration of the neu- 
tral hydroxyl radical, OH. For low rotational levels the 2 77 
electronic ground state of OH is sufficiently described by 
Hund's case (a). Spin-orbit coupling results in the 2 77 3 / 2 
state lying ~139 cm -1 below the 2 II 1 / 2 state. Because of 
the non-zero orbital angular momentum of the unpaired 
electron in OH each total angular momentum, J, state 
is A-type doubled, resulting in two closely spaced (~1.5 
GHz), opposite parity levels in the ro- vibrational ground 
state. The application of an electric field readily mixes 
these states and for the symmetric f state results in an 
increase in energy with electric field ^3] • The most polar- 
ized sub-level of this state, denoted as | 2 77 3 / 2 , J = 3/2, 
mj = 3/2, f ), is decelerated by our Stark decelerator. We 
detect the presence of OH molecules through the tech- 
nique of laser induced fluorescence. A pulsed excitation 
laser tuned to 282 nm promotes the molecules along the 
A 2 U 1/2 (v = 1, J = 3/2) ^X 2 n 3/2 (v = 0, J = 3/2) elec- 
tronic transition, excited molecules subsequently decay 
primarily (73%) back to the ground electronic state along 
the A 2 U 1 / 2 (v = 1) -^X 2 IJ 3 / 2 (v = 1) pathway at 313 nm. 
This fluorescence is collected and imaged onto a photo- 
multiplier tube. By varying the time of the laser pulse, we 
sample the OH molecules at a single location at different 
times to extract time-of-flight (ToF) information. Care- 
ful measurements of the total molecule numbers before 
and after the Stark decelerator, and thus determination 
of the overall decelerator efficiency as a function of de- 
celerator/supersonic beam operation parameters, enables 
information pertaining to the phase space matching of the 
source to the decelerator to be obtained. We have devel- 
oped a simple model along with results of Monte Carlo 
simulations that provide excellent agreement with the ob- 
served behavior. The remainder of the article is organized 
as follows. Section 2 briefly describes the operation of a 
Stark decelerator (for a more complete review of Stark 
deceleration see [T5UTS) ) as well as develops the simple 
intuitive model for phase space matching. Section 3 de- 
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Fig. 1. Schematic of the Stark decelerator, displaying the 
pulsed valve and discharge plates, the molecular beam skim- 
mer, the electric hexapole, and the electrode stages. The elec- 
trode stages alternate orientation (vertical - horizontal) as 
shown in the figure. The spatial locations indicated by arrows 
correspond to the locations where molecule number is mea- 
sured to determine decelerator efficiency. 



tails our experimental results and compares them to the 
model of Section 2 to provide confidence in using it for op- 
timization, while Section 4 is reserved for conclusions and 
suggestions for optimal operation of a Stark decelerator. 



2 Stark Deceleration 

Figure 1 shows a schematic of our Stark decelerator, which 
has been described elsewhere [El- Hydroxyl radical molecules 
produced by the discharge of H2O in Xenon (~1:99 mix- 
ture) undergo a supersonic expansion and are subsequently 
skimmed to separate the low vacuum source region (~10 -4 
torr) from the high vacuum region (~10 _T torr), which 
contains the decelerator electrodes and the applied high- 
voltage electric fields. After skimming, the molecules are 
focused by an electrostatic hexapole field to provide trans- 
verse coupling into the Stark decelerator. Once the molecules 
are coupled into the Stark decelerator, the slowing process 
begins. The Stark decelerator is constructed of 69 slowing 
stages spaced 5.475 mm apart with each stage comprised 
of two cylindrical electrodes of diameter 3.175mm sepa- 
rated axially by a 5.175 mm and oppositely biased at high 
voltage (±12.5 kV). Successive stages are oriented orthog- 
onally to each other to provide transverse guiding of the 
molecular beam. The geometry of the slowing stages pro- 
vides an electric field maximum between the electrodes 
with the field decreasing away from the electrode cen- 
ter as seen in Figure 2(a). Therefore, a weak- field seek- 
ing molecule travelling longitudinally down the decelera- 
tor will be decelerated as it moves into the region between 
the electrodes and will remain decelerated if the high volt- 
age is removed before the molecule passes through the 
electrode's field. In order to minimize the number of high 
voltage switches used, all the like-oriented electrodes of 
the same polarity are connected to one switch (4 switches 
total for the entire decelerator). 

As a consequence of this switch minimization, when 
the voltage at each slowing stage is grounded and the next 
stage is turned on, the electric field is not completely re- 
moved, but rather becomes the field generated by the al- 
ternate set of electrodes (i.e. the solid versus dashed curves 
in Figure 2(a)). Thus, the decrease in the molecular kinetic 
energy, ZAKE, is given as the difference in potential en- 
ergy generated by the two sets of electrodes and is shown 
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Fig. 2. Phase stable operation of the decelerator. (a) Longi- 
tudinal Stark energy potentials generated by the two sets of 
electrodes, where the solid line represents the potential from 
one set of electrodes and the dotted line from the other set. (b) 
Kinetic energy loss per stage (AKE) experienced by molecules 
from switching between the two potentials given above. Dashed 
line corresponds to a sine function approximation of the poten- 
tial energy loss. The shaded region corresponds to the phase 
stable area when decelerating at 45°. 



as a function of synchronous molecule phase angle, </>o, in 
Figure 2(b). The synchronous molecule phase angle, </>ois 
defined as: 

a 180 ° m 
0o = — j-z, (1) 

where L is the distance between two adjacent stages and 
the zero position of z is defined to be exactly between 
two adjacent stages. Thus, switching at 0o = 0° corre- 
sponds to no net change in molecular kinetic energy, while 
maximum deceleration (acceleration) occurs at 0o = 90° 
(-90°). Though synchronous molecule phases between </>o 
= 0° and 180° lead to deceleration, only 0° < O < 90° 
results in stable slowing because it is only on the positive 
sloping AKE curve that molecules ahead (behind) of the 
center of the bunch are decelerated more (less) than the 
average deceleration. To further investigate the dynam- 
ics of the molecules within this stable packet about the 
synchronous molecule (i.e. the center of the bunch) it is 



useful to define the excursion of a molecule from the syn- 
chronous molecule position as A(j> = <j) — 4>o- Then with 
the use of Newton's second law and the sine function fit 
to the change of kinetic energy as function of 0o (Figure 
2(b)) it is trivial to write |T7IT5| : 



d 2 A(\> t ttC 



Ma 



^(sin(Z\0 + 0o ) - sin(0o)) = 0, (2) 



rriL 2 



where CMax is the maximum kinetic energy change per 
slowing stage, t is the time coordinate, and m is the molec- 
ular mass. This equation is that of a harmonic oscillator 
with its equilibrium position offset (represented by 0o). 
Thus one expects non-synchronous molecules to oscillate 
around the synchronous molecule position inside an asym- 
metric oscillator potential. From numerical integrations of 
Eq. (2) and the first time derivative of Eq. (1) one can 
solve for the stable and unstable regions of phase space, 
shown as a function of the synchronous molecule phase 
angle in Figure 3 for OH under typical decelerator opera- 
tion. The most important feature of Figure 3 is the rapidly 
decreasing area of stable evolution (i.e. region bounded by 
the separatrix). This decrease in stable area is easily un- 
derstood by analogy to the pendulum driven by a constant 
torque. As the torque is increased (i.e. the synchronous 
phase angle is increased) the equilibrium position of the 
pendulum is pushed toward the apex, therefore the ampli- 
tude of oscillations that result in stable oscillatory motion 
is reduced. This behavior is responsible for the separatrix 
possessing only one cusp for non-zero 0q. It is clear that 
the number of molecules accepted by a decelerator operat- 
ing at a specific synchronous molecule phase angle is then 
given as an area integral inside the separatrix weighted 
by the supersonic beam distribution at the decelerator's 
entrance. 

This phase space "bucket" loading is illustrated graph- 
ically in Figure 4. In the first panel of this figure, ToF data 
taken at the decelerator entrance (location "a" indicated 
in Figure 1) are shown along with a hypothetical ideal 
supersonic beam of OH molecules (v^ center — 

300 m/s, 

v z = 30 m/s, where v ZyCente r is the pulse mean speed and 
Av z is the full width at half maximum of the distribu- 
tion (FWHM)). For this graph the supersonic molecular 
beam parameters were varied by tuning the discharge ini- 
tiation time as described previously ^BJ resulting in the 
ability to input two vastly different molecular beams into 
the decelerator. The first distribution utilized has v z center 
= 415 m/s with Av z of ~ 90 m/s, while the second dis- 
tribution is centered about 350 m/s with a spread of ~ 
80 m/s. For these two experimental distributions and the 
idealized case, Monte Carlo simulation results are used to 
construct graphs of the longitudinal phase space occupied 
by the molecules in the remaining panels of the figure. 
Each gray point in these graphs represents the location in 
phase space of a sample molecule as predicted by our simu- 
lations at the time of loading into the decelerator. The two 
dashed curves on each graph are projections of the phase 
space distribution onto the respective axes. Overlaid on 
these graphs are the separatrix curves for operation at </>o 
=0°, the so-called bunching condition, and for slowing the 
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Fig. 3. Phase stable area (separatrix area) versus phase angle 
for 0o = 0°, 30°, and 60°. As the phase angle is increased, 
the stable area is reduced; equivalently, the stable longitudinal 
velocity width narrows, and assuming the source completely 
populates the phase stable region, decreases the number of 
molecules in the stable packet. 



molecules to rest at the exit of the decelerator. From these 
graphs, the importance of creating a cold (i.e. small Ay) 
source of molecules as an input for a Stark decelerator is 
clear; with a smaller Ay the molecular pulse spreads less 
in route from generation at the source to the decelera- 
tor entrance, and thus is more efficiently coupled into the 
decelerator in both the spatial and velocity dimensions. 
The beneficial effect of lowering the pulse mean speed is 
also made evident; as the mean speed drops, the required 
phase angle for slowing is reduced and thus the separatrix 
area increases. 

As aforementioned, the effects of supersonic beam and 
Stark decelerator parameters on overall efficiency can be 
estimated from knowledge of the separatrix overlap with 
the supersonic beam phase space. Approximating the area 
enclosed by the separatrix as that of the corresponding 
rectangle whose height (width) correspond to the separa- 
trix velocity (spatial) bounds and ignoring the minimal 
convolution due to the short free flight of the molecules 
after creation, the overall longitudinal slower efficiency is 
given as: 



V Longitudinal 

XN[Z, Az, ZBoundsiM], 



(3) 



with 



N(x, Ax, X Bounds) 



/>XV pper -bound (4>o) 
X Lower bound(4>o) 



f(x, Ax)dx 



ST 00 f(x,Ax)dx 



(4) 



where z refers to the dimension along the decelerator's 
axis (longitudinal axis) and f(x, Ax) is the molecular dis- 
tribution in the appropriate dimension. The integration 
bounds correspond to the separatrix maxima along their 
appropriate dimension and are an explicit function of the 
synchronous molecule phase angle. From Figure 2(b) we 
see that for 0o = 45°, non-synchronous molecules with 
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Fig. 4. Longitudinal bucket loading. Panel (a) represents 
time-of-flight measurements taken at the decelerator input (po- 
sition "a" in figure 1) for different molecular beam operating 
conditions. Filled circles (filled squares) represent operating 
with a distribution with a mean speed of 415 m/s (350 m/s) 
and a full width at half maximum spread of 95 m/s (85 m/s), 
while the solid line is an idealized molecular beam of 300 m/s 
mean speed with 10% velocity spread. Panels (b), (c), and (d) 
represent the longitudinal phase space at the beginning of the 
Stark deceleration as derived from our Monte Carlo simula- 
tions. In panel (b) the dashed separatrix represents bunching 
for 415 m/s, while the larger solid line is the separatrix for 
bunching at 370 m/s. Inside the 370 m/s bunching separatrix 
is shown the separatrix for slowing a packet of molecules to 
25 m/s. In panel (c) the separatrix are shown for bunching at 
325 m/s and slowing to 25 m/s molecules from the 350 m/s 
distribution, while panel (d) shows bunching at 300 m/s and 
slowing to 25 m/s for the hypothetical molecular beam. 



phase angles between 45° and 135° will be slowed more 
than the synchronous molecule (shaded region). In gen- 
eral, we see the maximum stable forward excursion is then 



given as: 



Max,+ 



(<fo) = 180° - 2<h 



(5) 



Unfortunately, a closed form solution for the maximum 
backwards excursion, A(j)Max,-(4>o), does not exist, how- 
ever, it can easily be found numerically or estimated by 
the negative of Eq. (5) since both the maximum stable 
forward and backward excursions have the same average 
derivative as a function of 0r> From Eq. (2) we find the 
work done in bringing a molecule starting at Ac/) m ax, + (fa) 
with the same speed as the synchronous molecule (i.e. at 
the separatrix cusp) back to the synchronous molecule po- 
sition (i.e. A(j) = 0), W, as: 



W 



-Cmcl 



f 



(sm(Acj) + 0o ) — sm((j)o))dA<^ 

) 

(6) 
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Clearly, this molecule will possess the maximum stable 
velocity and we thus have the separatrix velocity bound 



Au Max (M = 2\ — — (cosOo) - - 0o)sin(0 o )). 

(7) 

Now, assuming the molecular distributions are ade- 
quately described by the form: 



Ax) = exp 



% ^center 



2VM2) 



(8) 



where Ax refers to FWHM and x cen ter to the value at 
which the distribution is centered we see Eq. (3) becomes: 



V Longitudinal — 
1 rz? n I (Vz,design ~ Vz, center) + Av Zi Max(<t>o) 



AErf 



n I ( v z, design ~ Vz, center) — Af z ^M ax (<t>o) H — 77^" li 

~ hr J I 2^7 v ln { 2 ) J\ 

( (Zdesign ~ ^center) + ^A(j)Max,+ (00 ) r, 

x W \ aT^ VH2) J 

t— i r f (^design ~ Z center ) — ^ A(/)Max r, /n x V 

~ Er f I y/ln(2) I] 

Where Erf refers to the Error function, and v z Resign 
and ^design have been introduced to account for the possi- 
bility of designing the decelerator pulse sequence to select 
molecules which are not at the peak of the molecular dis- 
tribution. If molecules are loaded at the peak of the dis- 
tribution (i.e. V z, design ~ Vz, center and ^design ~ ^center) 

and we approximate A^Max-iM = -^</Wax,+ (</>o), we 
recover the simple result of our earlier work [T5] : 



TjLongitudninal — Erf 



Av. 



z,Max{(P0) 



(io) 



xErf 



Ax 
2 



which is useful for estimating expected decelerator ef- 
ficiency. Plots of Eq. (8) with loading at the peaks of the 
distributions are shown in Figure 5 versus both Av z (Fig- 
ure 5(a)) and v z , center (Figure 5(b)), where the value of 
Ac/) Max, -((fro) nas been found numerically. In these graphs 
the slowing phase angle was chosen to bring the molecules 
to rest at the decelerator's exit. Again the importance of 
a low central velocity and small spread about this mean is 
evident. The most striking feature of these graphs is the 
flattening of the longitudinal slower efficiency for small 



velocity spreads. In fact for a sufficiently narrow velocity 
spread the mean speed of the pulse becomes to some de- 
gree "unimportant" as the supersonic beam emittance fits 
entirely inside the decelerator acceptance for most oper- 
ating conditions. It is important to note that similar ef- 
fects can be achieved by changing the characteristics of the 
Stark decelerator. By increasing the stage-to- stage spacing 
one enlarges the spatial bounds of the separatrix resulting 
in more decelerated molecules; however, increasing the dis- 
tance without increasing the applied electric field will re- 
sult in worse transverse guiding and a smaller stable phase 
space area transverse to the molecular beam axis. Increas- 
ing only the applied electric field results in the ability to 
decelerate molecules with larger velocity spreads relative 
to the synchronous molecule speed, and thus increases the 
overall slower efficiency. However, the technical challenges 
associated with going to voltages higher than the present 
operating conditions are large. Therefore a more practi- 
cal alternative maybe to increase the number of slowing 
stages to reduce the required phase angle. 

Efficient operation of a given decelerator /supersonic 
beam combination occurs when Eq. (8) is maximized at 
the desired final speed. Since the molecular distribution 
and the desired final speed are assumed fixed the only pa- 
rameter which can be varied is v z ,desigm the initial speed 
for which the decelerator pulse sequence is designed. By 
decreasing v z , design the required slowing phase angle is 
decreased and the separatrix area increased. However, if 
v z, design is lowered too far below v ZjC enter the molecular 
beam phase space overlaps poorly with the decelerator in- 
put separatrix and slower efficiency suffers. Thus, the most 
efficient slower operation will occur for a v z , design below 
Vz, center by an amount that depends on the Av z of the 
pulse and the desired final speed (typically 1 to 10 per- 
cent of v z , center for slowing to rest). Figure 6 shows the 
results of maximizing Eq. (8) for the given experimental 
distributions as well as the idealized source. In this graph 
the amount of velocity detuning from v Zi center for most 
efficient operation is shown as a function of desired final 
speed. For final packet speeds close to the initial pulse 
speed (right part of graph) it is most efficient to utilize a 
phase angle of zero (bunching) because of its large accep- 
tance and vary v z Resign - However, as the final velocity is 
decreased, this procedure eventually leads to a mismatch 
between the supersonic beam emittance and the decelera- 
tor acceptance, and it becomes more efficient to decelerate 
molecules nearer the center of the supersonic beam distri- 
bution, leading to the sharp "transition" peaks. The posi- 
tion of this transition for our decelerator is at a detuning 
equal to approximately 15% of the pulse spread and is evi- 
dent in the figure by the sudden drop in velocity detuning. 
As slower molecules are desired, the effect of the decreas- 
ing acceptance of the decelerator (with increasing </>o) is 
combated by selecting slower initial molecules, which re- 
quire less slowing. This effect is evident by the increase of 
detuning with slower final speeds. The gain in efficiency 
by optimizing decelerator operation according to Eq. (8) 
depends sensitively on both the molecular beam and de- 
celerator parameters, but in the case of our experimental 
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Fig. 5. Simple model for longitudinal bucket loading. Panel 
(a) shows the effect of varying the supersonic beam's velocity 
spread on longitudinal efficiency for several initial speeds, while 
panel (b) shows the longitudinal efficiency as a function of 
initial pulse speed for several velocity spreads. 



conditions ranges from only a few percent to well over a 
factor of 100 improvement [T9] . 

While Eq. (8) describes the efficiency of decelerator 
operation in the longitudinal dimension, it says nothing 
about evolution transverse to the decelerator axis. In this 
dimension the molecular distribution is centered about 
zero in both space and velocity thus Eq. (3) becomes: 
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Fig. 6. Decelerator optimization condition. The optimum de- 
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sonic beam distributions. 
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^Transverse 

( Ax{cj> ) 



Erf ( %^v^ 



I Ax 

\ 2 



VM2) )Erf 



Av y ,Max(4>o) 



*Er f ( 



\ Ay 

\ 2 



(12) 



The spatial acceptance of the decelerator is set by the 
decelerator rod spacing, thus Ax and Ay ~ 1 mm. Be- 
cause the transverse evolution of molecules through the 
Stark decelerator is a complicated dynamical process that 
depends sensitively on the operating conditions it is not 
possible to give a simple expression for Av Xi Max(<f> ) an d 
Av Xi Max((/>o)i however, for most deceleration experiments 
to date it can be estimated as ^3 m/s for most slower 
conditions with minimal error 16 , 20 . The behavior of Eq. 
(11) is analogous to that of Eq. (8), as the velocity (spa- 
tial) spread of the source is lowered or the maximum sta- 
ble velocity (position) is increased by a decelerator change, 
the efficiency of the slower grows. It is the aim of proper 
hexapole focusing to "mode-match" the supersonic beam 
into this transverse acceptance set by the decelerator. 



^Transverse 

= N[0, Av x , V x ,Bounds(M}N[0, Ax, X Bounds(4>o)] 
N[0, Av y , Vy,Bounds((l)o)}N[0, Ay, y Bounds^)], (H) 



where x and y refer to the dimensions perpendicular to 
the decelerator's axis. Assuming distributions of the form 
of Eq. (7) we have: 



3 Experiment 

Time of flight measurements of OH molecules taken at 
the exit of the decelerator (position "b" in Figure 1) are 
shown in Figure 7 for various decelerator operating condi- 
tions. Open circles represent operation in transverse guid- 
ance mode, where the decelerator array is biased at high 
voltage, but un-switched resulting in a potential minimum 
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Fig. 7. Molecular packets at the decelerator exit (position "b" 
in Figure 1) for varying phase angles. Peaks arriving later in 
time are the result of slowed molecules. From comparison of 
the number of molecules in these packets to the total inputted 
molecular number, decelerator efficiency can be determined. 



along the decelerator axis. This mode of operation does 
not affect the molecules longitudinal speed, but provides 
transverse guidance through the length of the decelerator. 
The remaining traces show the stable packet of molecules 
produced when switching the decelerator with non-zero 
phase angle. As the slowing angle is increased the packets 
arrive later in time signifying their smaller mean speed. 
The spreading of the ToF pulse due to free flight and 
convolution with the detection window is also evident. 
Integration of the deconvolved ToF pulse reveals the to- 
tal molecular number in each phase stable packet, and is 
shown in Figure 8, where the data points have been nor- 
malized to the bunching condition. In this graph the total 
molecular number is plotted versus the final speed of the 
molecular bunch for the three different operating condi- 
tions represented in Figure 4. The filled circles represent 
data taken for slowing sequences designed for an initial 
speed of 415 m/s. For this data the supersonic beam was 
operated to yield a distribution centered at 415 m/s with 
a 90 m/s FWHM. The filled squares in Figure 8 represent 
loading molecules from the same distribution into a phase 
bucket designed for 370 m/s initial speed. For the data 
represented by the filled triangles, the supersonic beam 
was operated to yield a pulse centered at 350 m/s with an 
80 m/s spread. In this data, the decelerator was operated 
for a 325 m/s initial speed. The solid lines in this graph 
represent the simple theory of Section 2, while the open 
data points correspond to the results of our Monte Carlo 
simulations. Despite the assumptions made in deriving Eq. 
(8) we see the simple theory of Section 2 is quite accurate. 
The only noticeable discrepancy occurs for the 415 m/s 
curve, where it appears that deviation in the experimen- 
tal bunching point (which is used to normalize the graph) 
is responsible for the offset. From this graph we see that 
operation at 325 m/s is clearly more efficient at produc- 
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Fig. 8. The stable molecule number normalized to (f)o— 0° 
for three cases of operation. Circles represent operation with a 
supersonic beam centered at 415 m/s with a 90 m/s FWHM 
and the deceleration sequence designed for molecules at 415 
m/s. Squares represent a deceleration sequence designed for 
370 m/s molecules from the same beam distribution. Triangles 
are for deceleration designed for 325 m/s molecules from a 
distribution centered at 350 m/s with an 80 m/s spread. Filled 
points represent data, while open points represent results from 
the Monte Carlo simulation. The solid lines are the results of 
the simple model of Section 2. 



ing slow molecules. Again, the importance of developing 
a slow cold source as the input for efficient Stark deceler- 
ation is clear. 

By comparing the total molecule number observed in 
each packet to the total number of molecules inputted into 
the decelerator (Figure 4 (a)) the overall efficiency of the 
Stark decelerator can be determined and is shown in Fig- 
ure 9 for the 370 m/s case as filled squares. Overlaid on 
this graph is the simple model of Section 2. For this the- 
ory the transverse parameters Z\x and Ay = 1.5 mm and 
Av x and Av y = 12 m/s were used in Eq. (11) |STj . It is 
remarkable to note that this simple theory predicts the 
overall slowing efficiency quite accurately with no free pa- 
rameters. Clearly, this level of agreement with our simple 
model gives confidence in using it as an optimization tool. 



4 Conclusions 

We have developed a simple model of the process of phase 
space matching between the supersonic beam emittance 
and the decelerator acceptance. The ability of this model 
to predict with accuracy the efficiency of the Stark de- 
celeration process gives confidence in its use as an opti- 
mization tool. While optimization according to this model 
gives small gain for supersonic sources with low mean ki- 
netic energy and small velocity spreads, it can lead to dra- 
matic improvement for molecular beams that have either 
a high kinetic energy, a large velocity spread, or both. As 
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Fig. 9. Stark deceleration efficiency. Efficiency of the Stark 
decelerator as a function of final speed for slowing designed 
for an initial speed of 370 m/s taken from the distribution 
centered at 415 m/s with a 90 m/s FWHM. Filled squares are 
data points obtained from the results of Figure 8 normalized 
by the results of trace (a) in Figure 4. The solid line is the 
model of Section 2 including both longitudinal and transverse 
efficiencies. 



the technique of Stark deceleration is extended to more ex- 
otic molecular species, which because of their large mass, 
internal modes or process of creation result in non-ideal 
supersonic beams for input into a Stark decelerator, max- 
imization of the decelerator operation efficiency will be 
essential. 

Clearly, it is the aim of the experimenter to maximize 
the product of Eqs. (8) and (11) for the most efficient de- 
celerator operation. The efficiency of the decelerator de- 
pends sensitively on both the supersonic beam and Stark 
decelerator parameters. The importance of developing a 
cold source as the input for any Stark deceleration ex- 
periment cannot be overstated since this is the only real 
cooling in a slowing experiment and, for a given deceler- 
ator, sets the attainable efficiency. While increasing the 
dimensions of slowing stages and their separations to al- 
low deceleration of more molecules is a viable way to in- 
crease the deceleration efficiency, it is met with increas- 
ing technical challenge as the electrode voltage must also 
be scaled. Perhaps a more reasonable alternative is to in- 
crease the number of slowing stages, thus dropping the 
amount of energy required to be removed at each stage, 
resulting in a larger phase space acceptance. A reasonable 
design goal is to include enough stages in the decelerator 
so that molecules can be decelerated to the desired final 
speed with a phase angle of 0o < 45°. Also, positioning 
the source as close as possible to the input of the decelera- 
tor is of importance to provide efficient spatial coupling in 
the longitudinal dimension. Finally, for a given molecular 
beam/decelerator combination the optimum phase angle 



and initial pulse design speed, v z center , for slowing can 
be found from Eq. (8). 
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